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Abstract. We have observed the continuum emission of the interacting galaxy NGC 4631 at AA 870 /im and 
1.23 mm using the Heinrich-Hertz-Telescope on Mt. Graham and the IRAM 30-m telescope on Pico Veleta. We 
have obtained fully sampled maps which cover the optical emission out to a radius of about 7' at both wavelengths. 
For a detailed analysis, we carefully subtracted the line contributions and synchrotron and free-free emission from 
the data, which added up to 6% at 1.23 mm and 10% at 0.87 mm. We combined the flux densities with FIR 
data to obtain dust spectra and calculate dust temperatures, absorption cross sections, and masses. Assuming a 
"standard" dust model, which consists of two populations of big grains at moderate and warm temperatures, we 
obtained temperatures of 18 K and 50 K for the both components. However, such a model suffers from an excess of 
the radiation at A 1.23 mm, and the dust absorption cross section seems to be enhanced by a factor 3 compared to 
previous results and theoretical expectations. At large galactocentric radii, where the galaxy shows disturbances 
as a result of gravitational interaction, this effect seems to be even stronger. Some possibilities to resolve these 
problems are discussed. The data could be explained by a very cold dust component at a temperature of 4-6 K, 
an increased abundance of very small grains, or a component of grains with unusual optical properties. We favour 
the latter possibility, since the first two lead to inconsistencies. 
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, !_J 1. Introduction 



The efficiency and time-scale of star-formation processes 
in galaxies depend strongly on the amount, distribution, 
and composition of the interstellar matter. One compo- 
nent of the ISM whose properties are still not well known 
is interstellar dust. 

Since most of the dust in galaxies is cold, with tem- 
peratures of < 20 K, infrared observatories like IRAS are 
almost blind for the gross amount of dust, and can only 
observe the warm dust component. The major fraction of 
the interstellar dust radiates mainly in the sub-mm range, 
at wavelengths > 100 /im. This emission was hardly ac- 
cessible in the past, due to the lack of good sub-mm tele- 
scopes at sites with sufficiently good atmospheric condi- 
tions. This situation has improved in the last few years: 
Instruments like the bolometer arrays installed at the 30- 
m telescope on Pico Veleta or - more recently - SCUBA 
at the JCMT and the 19-channel array at the Heinrich- 
Hertz-Telescope have produced important results for the 
investigation of this cold dust component. Also, the ISO 



satellite could routinely measure the peak of the cold dust 
emission from grains with temperatures ranging down to 
10 K (e.g. Popescu et al.EEEl- 
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However, dust properties in external galaxies are still 
not well determined because of gaps in the spectral cover- 
age. Another difficulty in the analysis and interpretation 
of continuum data is the possible contribution of molecu- 
lar lines to the total flux measured with the (broad-band) 
bolometer arrays, especially in the usually used atmo- 
spheric windows around AA870/im and 1.23 mm which 
contain the strong CO(3-2) and (2-1) lines. This situation 
currently improves as more and more objects are mapped 
even in the higher CO transitions. 

This current paper is the first in a series that reports 
on observations of the thermal emission of cold dust in 
selected galaxies. NGC 4631 is a nearby (D — 7.5 Mpc, 
e.g. Golla & Wielebinski [T9M| galaxy which is close to 
edge-on. While the position angle of the large-scale (cm-) 
radio continuum emission is about 86°, that of the inner 
disk is closer to 82°. A value of 84°, which we assume 
throughout this paper, seems to be most appropriate for 
the mm/sub-mm emission out to a radius of a few arcmin- 
utes. NGC 4631 is embedded in a small gravitationally in- 
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Table 1. Some basic parameters of NGC 4631 as obtained 
from the literature. 



Type 


Sd 


(de Vaucouleurs et al. 1991 ) 


Position: 




(Young et al. H^5l 


R A [2000] 


12 h 42 m 07!65 




Dec[2000] 


32°32'27!9 




Distance 


7.5 Mpc 


(Golla & WielebinskirKMt 


Pos. Angle 


84° 




Inclination 


86° 




M (H i) 


5 x 10 9 M Q 


(Rand[HMl 


M(H 2 ) 


1.2 x 10 9 Mq 


(Golla & WielebinskinMll 


IRAS flux densities: 


(Young et al. HMfll 


12 jum 


5.1 ± 1.1 




25 /im 


8.8 ± 1.8 




60 /im 


90.7 ± 18.2 




100 /im 


170.4 ± 34.1 





teracting group with two neighboring galaxies. The dwarf 
elliptical galaxy NGC 4627 is located only 3' northwest 
of the nucleus of NGC 4631, and another edge-on spiral, 
NGC 4656, is located about 30' to the southeast. This in- 
teraction has created several prominent H I bridges and 
spurs (Rand 1994). Since galaxy interactions play an im- 
portant role in the evolution of galaxy systems and the 
ISM content in galaxies, the investigation of this object 
can provide useful information about many processes in 
the ISM. 

NGC 4631, which was classified as "mild starburst" by 
Golla & Wielebinski (1994J, may be in a late stage of its 
interaction, where the central star formation (triggered 
by molecular inflow due to the gravitational forces) has 
already ceased and an energetic outflow as well as a huge 
radio halo have been formed (Ekers & Sancisi 1977; Rand 

mm . 

The observations of the dust component of NGC 4631 
started with the IRAS satellite (Rice et al. 119881 Young 
et al. IT§5§|) . While the IRAS observations were sensi- 
tive mainly to the warm dust present in this galaxy, the 
A 1.3 mm map of Braine et al. ( 1995 ) proved the existence 
of a significant amount of cold dust in the central area. 
The inner 2'.h were also observed at A 850 /im by Alton et 
al. I|1999f) . Neininger & Dumke |T999| presented a more 
extended map at A 1.23 mm (which was obtained from the 
same data as shown in this paper) and detected inter- 
galactic cold dust which was pulled out of the disk by 
the gravitational interaction of NGC 4631 with its neigh- 
bours. Here we present a more extended map at A 870 /im, 
covering the disk out to a radius of 7'. 

Some basic parameters of NGC 4631 are compiled in 
Table □ 

2. Observations and data reduction 

We observed the galaxy NGC 4631 at wavelengths 
of AA 870 /im and 1.23 mm, using sensitive 19-channcl 
bolometer arrays developed by E. Kreysa and collabo- 
rators at the Max-Planck-Institut fur Radioastronomie, 
Bonn. 



2.1. HHT observations 

The observations at A 870 /im were carried out at the 
Heinrich-Hertz- Telescope 1 (Baars et al. I1999|l . located on 
Mt. Graham, Arizona, during three observing sessions be- 
tween March 2000 and January 2002, using a 19-channel 
bolometer array installed as facility instrument. The 19 
channels of this bolometer are located in the centre and 
on the sides of two concentric regular hexagons, with an 
apparent spacing between two adjacent channels (beams) 
of 50". The central frequency of the bolometer is about 
345 GHz (the highest sensitivity is reached at 340 GHz) , 
and the instrument is sensitive mainly between 310 and 
380 GHz. 

In order to calculate the atmospheric zenith opacity, 
we made skydip observations every 45 to 90 minutes, de- 
pending on the atmospheric stability. These yielded at- 
mospheric opacities at our observing frequency which var- 
ied between 0.3 and 0.9 for the several observing sessions. 
For calibration purposes we have also performed mapping 
and on-off measurements of various planets (mainly Mars 
and Saturn) during the observations. These measure- 
ments yielded a conversion factor from observed counts to 
mjy/beam of 0.8 - 1.1 mJy beam" 1 count -1 . This vary- 
ing conversion factor is due to (and partly corrects for) 
varying atmospheric condition at the day of the observa- 
tions and the uncertainties in the opacity calculation. The 
beamwidth at this frequency is ~ 23". 

All maps were observed in the Az-El coordinate sys- 
tem, by scanning along Azimuth and with data acquisi- 
tion every 0.5 s. During the observations, the subreflector 
was wobbled at 2 Hz in azimuth, with a beam throw be- 
tween 60" and 200". If necessary, the starting point of 
each subscan was shifted by a few arcseconds with re- 
spect to the preceding one, in order to ideally place each 
individual coverage relative to the target source. This as 
well as different map sizes (typically around 450" x 300") 
ensured that each subscan covers the galaxy and a suf- 
ficient amount of blank sky on either side, in order to 
facilitate baseline subtraction and a correct restoration of 
the double-beam maps. 

2.2. Pico Veleta observations 

The observations at A 1.23 mm were carried out in March 
1997 at the IRAM 30-m telescope on Pico Veleta, Spain. 
The central frequency of the used bolometer array is 
243 GHz, the bandwidth about 70 GHz. The array layout 
is similar to the one at the HHT, with individual chan- 
nels separated by 20". The beamwidth at the observing 
frequency is 11". 

The sky opacity at the observing frequency was moni- 
tored with sky dips and varied between 0.1 and 0.3 dur- 
ing the five days of the observing run, but was stable 
within ±0.02 for each individual day of the observations. 

1 The HHT is operated by the Submillimeter Telescope 
Observatory on behalf of Steward Observatory and the MPI 
fur Radioastronomie. 
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Fig. 1. Continuum map of NGC 4631 at A 0.87 mm, overlaid on an optical image extracted from the Digitized Sky- 
Survey. The map is smoothed to an angular resolution of 30"; the beam size is indicated by the filled grey circle in 
the lower left corner. The rms noise in the map is about 20mJy/beam area, and contour levels are 40, 80, 120, 160, 
200, 250, 300mJy/beam area. The dark spot in the optical (greyscale) image 3' northwest of the central area is the 
neighbouring galaxy NGC 4627. 



We mapped Mars every night to determine the absolute 
flux density scale; from these measurements we obtained 
a conversion factor of 0.24 - 0.28 mJy beam -1 count -1 . 

The observing procedures were similar to the HHT 
observations. The maps were observed in the Az-El co- 
ordinate system, but with a lower scanning velocity and 
smaller subscan separation in order to ensure full sam- 
pling. The subreflector was wobbled at 2 Hz with a beam 
throw of 45". 

2.3. Data reduction 

The data reduction for the HHT and the 30-m data was 
performed with the NIC program of the GILDAS software 
package. In addition, we used the MOPSI program to con- 
firm (and adjust, if necessary) the zenith opacities calcu- 
lated by NIC from our skydip measurements, and also to 
estimate flux densities for the observed planets in order 
to calibrate the data in astronomical units. After baseline 
subtraction and the elimination of spikes in each single 
coverage, the atmospheric noise, which is highly correlated 
between the individual channels (especially at 870 /mi), 
was subtracted. The maps were gridded, restored (using 
the EKH algorithm), converted into the RA-Dec system 
and finally combined (with an appropriate weighting) to a 
single map for each of the two wavelengths. The zero-levels 
of the resulting maps were checked and carefully adjusted. 
Map features smaller than the telescope beam appearing 



in the final maps were filtered out using a Fourier filter 
technique. For further analysis, both maps were smoothed 
to a FWHM of 24". This was done in order to enable a di- 
rect comparison with the existing CO data (see Sect. 13. 3|) 
and other ISM components. The rms noise level in the fi- 
nal maps is 20mJy/beam at A 0.87mm, and 3mJy/beam 
at A 1.23 mm. 

3. Results 

3.1. Total flux densities 

To measure the total flux densities, we used a ring integra- 
tion method which we applied to the disk of NGC 4631. 
While the parameter choice for this ring integration had 
some influence on the results, its effect was much smaller 
than the uncertainties coming from the observations and 
the data reduction process, including the calibration. Also 
the actual rms noise in the final maps has a rather small 
effect. We estimate the error of the flux densities to be 
10% at A 1.23mm, and 15% at AO. 87mm. 

The total flux density of NGC 4631, including all disk 
emission, is So.87mm = 3.78 ± 0.57 Jy for the observing 
wavelengths of 0.87 mm. At A 1.23 mm, the total flux den- 
sity is more difficult to determine. The reason for this 
are the bridges and spurs which were detected at this 
wavelength, as well as the fact that the edge of the de- 
tected emission is more due to the edge of the mapped 
region rather than due to the edge of the dust distribu- 
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Fig. 2. Continuum map of NGC 4631 at A 1.23 mm, overlaid on an optical image extracted from the Digitized Sky 
Survey. The map is smoothed to an angular resolution of 20"; the beam size is indicated by the filled grey circle in the 
lower left corner. The rms noise in the map is about 3mJy/beam area, and contour levels are 10, 20, . . . , 80mJy/beam 
area. 



tion. Limiting the integration of the observed emission 
to the main disk of NGC 4631, we estimate Si. 23mm = 
2.18 ± 0.22 Jy. Note that we got these values from inte- 
grating over the same area and using maps at the same 
angular resolution (24") for both wavelengths. 

These values cannot be easily compared to previous 
flux density estimates, since our maps are the first to cover 
(almost) all of the galaxy's emission. Braine et al. (1995) 
observed the central 4' x 3' of NGC 4631 at a wavelength 
of 1.3 mm and obtained a flux density of 0.64 Jy. When 
we restrict our integration to the same area, we measure 
~ 0.9 Jy. This difference is only partly due to the slightly 
different wavelength. Another reason for their lower value 
might be the lack of sensitivity for the off-plane emission 
and therefore an imperfect baseline subtraction - note the 
negative areas north and south of the major axis and the 
overall smaller intensity values in Fig. 1 of Braine et al. 

Bendo et al. (2005|) used SCUBA archive data to es- 
timate the flux density at A850/im and obtained 1.89 ± 
0.19 Jy and 0.54 ± 0.05 Jy for the central 135" and 45", 
respectively. For the same areas, we get values of 1.8 Jy 
and 0.5 Jy, which is in good agreement with those values. 

3.2. Morphology of the continuum maps 

The resulting maps of NGC 4631 at AA 0.87 mm and 
1.23 mm are shown in Figs. ^ and [21 overlaid on an op- 
tical image extracted from the Digitized Sky Survey. 



Similar to other ISM components, especially the CO 
molecular line emission, both maps show a double-peaked 
appearance in the central area, which points to a ring-like 
distribution with a diameter of about 1', or 2.2 kpc at the 
assumed distance of 7.5 Mpc. 

At larger radii, the emission is more patchy and fol- 
lows the major axis of the galaxy. The extent of the emis- 
sion in the disk is difficult to determine, since our map 
does not cover radii larger than 6' — 7' with a sufficient 
signal-to-noise ratio. At A 1.23 mm we find a significant 
amount of emission at distances of a few arcminutes away 
from the major axis, e.g. at (Ace, AS) = (250", -150") 
or (0", 120"), and thus several kpc above the plane. This 
emission is correlated with the Hi spurs found by Rand 
(1994J and is discussed in detail by Neininger & Dumke 
(1999). It is most likely interstellar material pulled out 
of the plane by the gravitational interaction of NGC 4631 
with its neighbouring galaxies. The total flux density of 
this extraplanar gas - as far as it is within our map edges 
- is about 0.5 Jy. This is one quarter of the flux density 
within the plane of NGC 4631, which is 2.18 mJy. 

At A 0.87mm we are unable to detect this intergalactic 
cold dust with sufficient significance. Only one area, lo- 
cated at (Aa, AS) = (60", 170"), is clearly detected. This 
feature is located on the Hi spur pointing towards the 
northeast as seen by e.g. Weliachew et al. (1978). None 
of the other extraplanar features of the 1.23 mm map can 
be seen above the noise in the sub- mm map. However, an 
integration of the map shows that a flux density of ~ 1 Jy 
originates in the off-plane features which can be identi- 
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fied in the 1.23 mm map, of which 400 mJy are due to the 
emission region described above, located northeast of the 
nucleus. The total flux density confined to the plane of 
NGC4631 and located within our map is 3.78 ± 0.57 Jy, 
so again the off-plane features contain about 25 % of the 
disk emission. 

The radiation measured by the bolometers contain 
some contributions which are not due to thermal dust 
emission. We calculate these non-dust contributions in the 
next subsection, before we continue with the discussion of 
the major axis distribution of important ISM components 
and the dust spectrum. 

3.3. Origin of the observed emission 

The broad band emission of the galaxies, measured with 
the bolometers, consists of four main components: thermal 
dust emission, free-free radiation from thermal electrons, 
synchrotron radiation from relativistic electrons, and the 
CO as well as some weaker lines which fall into the band- 
pass. As we are mainly interested in the thermal dust emis- 
sion alone, we have to determine the contribution of the 
other processes and to subtract them from the data. 

Since CO mapping observations of NGC 4631 exist, we 
can calculate -BcOj the contribution of the CO lines to the 
surface brightness measured with the bolometers, through 



Bco 



A^boi 



Jco 



(1) 



where the line intensity Jco = J T m b dv denotes the 
velocity-integrated main-beam brightness temperatures. 
From the integrated CO line flux Fco, the contribution 
to the total flux densities can be obtained via 



Sco — 



Ai/boi 



F, 



CO 



(2) 



Here the equivalent bandwidth of the bolometer can be 
calculated from the instrumental bandpass (Kreysa, priv. 
comm.) and is Afboi ~ 50 GHz for the CO(3-2) line and 
A^boi ~ 70 GHz for the CO(2-l) line. 

The contribution of the CO(3-2) to the 870 /im flux 
can be taken from data published by Dumke et al. (2001), 
and the CO(2-l) contribution to the 1.23 mm flux can be 
estimated from the data published by Golla & Wielebinski 



(1994). We should note that at both wavelengths the CO 
emission was observed with the same telescope as the cor- 
responding continuum emission, thus eliminating possible 
uncertainties due to beam characteristics. The total line 
contribution, including also 13 CO and other lines, may be 
somewhat higher. Based on typical line ratios in exter- 
nal galaxies and the sensitivity variation of the bolometer 
over the bandpass, we estimate that other lines add up 
to about 5% of the flux calculated from the 12 CO(3-2) 
line, and about 10% of the flux from the 12 CO(2-l) line. 
Taking this into account and entering all quantities in Eq. 
0J, it simplifies to 



Table 2. Flux densities of NGC 4631 at AA 0.87mm and 
1.23 mm. The non-dust contributions at both wavelengths 
are calculated as described in Sect. 13.31 



A 


5'obs [Jy] 


5i in e [mJy] 


Sff+sync [ m Jy] 


Sdust [Jy] 


0.87 mm 


3.78 ±0.57 


330 ± 30 


40 ± 10 


3.41 ±0.58 


1.23 mm 


2.18 ±0.22 


90 ±10 


50 ± 10 


2.04 ±0.23 



and 



B 



1.23mm 
line 



mJy] = 0.303 1, 



CO(2-l) 



[K km s 



(4) 



B° n 8 e 7mm [mJy] = 1.36 J CO (3-2) [Kkms" 1 ] 



(3) 



for the AA 870 /im and 1.23 mm emission, respectively, and 
a 24" beam. To correct for the line contribution in the 
continuum maps, we scaled the CO intensity maps ac- 
cordingly and subtracted these from the bolometer maps. 
For a comparison of the dust distribution at both wave- 
lengths with other ISM components (see next subsection) 
we smoothed all maps to a final resolution of 24" HPBW. 

For the line contributions to the total flux denity, we 
get S' 1 ° n 8 c 7mm = 330 ± 30 mJy and <S£ n 2 e 3mm = 90 ± 10 mJy. 
The fraction of this contribution to the observed flux den- 
sity is thus 9 % and 4 % for the two wavelengths of 0.87 mm 
and 1.23 mm, respectively. 

The fraction of the free-free and the synchrotron emis- 
sion can be estimated from radio continuum data at lower 
frequencies. Niklas et al. <|1995l I1997f) separated these two 
components in the radio spectra of a large sample of galax- 
ies. At a frequency of 10 GHz they found a thermal frac- 
tion of the emission of fth = 0.13 and a nonthermal spec- 
tral index of a n th = 0.78. With these values, we esti- 
mate a contribution of the free-free and synchrotron emis- 
sion at both wavelengths of 40 ± lOmJy (0.87mm) and 
50 ± 10 mJy (1.23 mm), which is of the order of 1 - 2 % of 
the total emission. Although the thermal fraction of the 
cm-emission might be underestimated because of galactic 
wind and cosmic ray propagation effects on the integrated 
radio spectrum (Werner 1988), the contribution of free- 
free and synchrotron radiation is still less than 80 mJy at 
both wavelengths even in the unlikely case that more than 
half of the emission in the cm-range is free-free emission. 
Therefore it is much smaller than the line contribution, 
especially at A 0.87 mm. Thus any uncertainties in the as- 
sumed thermal (free-free) fraction of the cm-emission and 
the non-thermal spectral index are negligible for the fur- 
ther data analysis, considering the total uncertainty of the 
flux density values of 10 - 15 %. 

The contributions of the line emission and the non-dust 
continuum radiation to the total flux densities at 0.87 mm 
and 1.23 mm are listed in Table 2. At this point we should 
note that the main uncertainties in the results are due to 
the rms noise in and the absolute calibration of the con- 
tinuum maps. The subtraction (or non-subtraction) of the 
non-dust contribution does not change our results quali- 
tatively, and the uncertainties in these contributions are 
negligible. 
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x offset [arcsec] 

Fig. 3. Distribution of various ISM components along the major axis of NGC4631. The thermal dust emission at 
870 /im and 1.23 mm, as obtained from the maps in Figs. ^ and [2] and corrected for line, synchrotron and free-free 
contributions, are shown as filled circles and white squares. The errors of both data sets are indicated in the upper 
left corner, while the long-dashed line gives their zero level. Also shown for comparison are the Hi profile (solid line), 
the CO (1-0) distribution (dashed line), and the A 6 cm radio continuum distribution (dotted line). The scale for these 
three data sets is shown at the y-axis on the right. 



3.4. ISM distribution along the major axis 

We determined the distribution of various components of 
the ISM along the major axis of NGC 4631, including the 
current results from AA0.87 and 1.23 mm (after subtrac- 
tion of the non-dust contributions), the Hi (Rand 1994), 
the CO(l-O) (GoUa & Wielebinski 11531 . and the A 6 cm 
radio continuum emission (Krause et al., in prep.). The 
result is shown in Fig. with an angular resolution of 
24" for all data sets. Here we averaged the data over three 
pixels (pixel size is 8") perpendicular to the major axis for 
each data point in order to increase the signal-to-noise ra- 
tio, and to account for emission which is close to, but not 
exactly on the major axis. The A 6 cm data are included as 
a tracer for the FIR emission, since at this wavelength we 
see mainly non-thermal radiation, for which the radio-FIR 
correlation is strongest. 

Similar to other edge-on galaxies, the dust emission 
resembles the distribution of molecular gas in the inner 
part of the galaxy. However, at radii where no molecular 
gas can be detected (r > 4'), there is still a significant 
amount of thermal dust emission, which in this case fol- 
lows the distribution of the H I. 

By comparing the distributions for the various wave- 
lengths of radio continuum emission, we note some signif- 
icant differences between the data sets. The central maxi- 
mum is only visible in the A 6 cm data, but not at mm/sub- 
mm wavelengths. Since the existence of an active nucleus, 
capable to produce excess cm emission, is rather unlikely 
in NGC 4631, this points to a steeper FIR-to-mm spec- 
trum and thus to higher dust temperatures in the very 



centre. Also the intensity of the "dust ring" is different 
at the eastern and western maximum, suggesting differ- 
ent dust properties in these two maxima. All distributions 
show a steep decrease in intensity at the eastern edge of 
the molecular ring, at a radius of 40" — 70", while the 
intensity decrease at the western edge is much shallower. 

At the molecular ring itself, at a radius of about 30", 
all continuum distributions show local maxima. Such a 
correlation between molecular gas and dust is expected 
because many molecules are formed on the surface of 
dust grains, and the dust shields the gas from starlight 
and prevents the molecules from being photodissociated. 
However, the H I distribution shows local minima at these 
positions. Here the interstellar hydrogen was probably 
transformed into molecules because of a local increase of 
ambient density. This is a direct consequence of gas flows 
due to a non-axisymmetric potential, which results from 
the gravitational interaction with the neighbouring galax- 
ies NGC 4656 and NGC 4627. 

While the CO intensities have dropped to relatively 
low values beyond the molecular ring, the Hi emission 
stays at significant levels. Also the tracers of interstellar 
dust (AA 1.23mm, 0.87mm, and 6cm, i.e. FIR) are still 
detected at larger radii. However, different correlations 
between the various ISM components can be observed. 
In the western disk, 3' — 4' away from the centre, the H I 
distribution shows a plateau-like feature. Here also the 
0.87 mm and 6 cm data show maxima in their distribu- 
tion, while we see only a small increase in the 1.23 mm 
data. In the eastern disk, on the other hand, the 1.23 mm 
and 6 cm distributions drop beyond a radius of 3', while 
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the 0.87 mm emission stays at a higher level, similar to 
the Hi. These significant differences in the distribution of 
the various dust tracers point to changing dust properties 
along the major axis of NGC 4631. This will be discussed 
m more detail in Sect. IP1 

4. Dust properties in NGC 4631 

4.1. The observed dust spectrum 

Table [21 lists the observed flux densities for both wave- 
lengths. In addition, it gives the non-dust contributions 
to the flux densities and the resulting values which are 
due to the thermal emission of dust alone. As discussed in 
Sect. 13.31 the line contribution of the measured continuum 
flux densities is 9% at AO. 87 mm, and 4% at A 1.23mm. 
The latter value is significantly below this value measured 
in other galaxies, either normal (e.g. NGC 5907, Dumke 
et al. lT$97|) . active (e.g. NGC 3079, Braine et al. lTM7|> . or 
starbursting (e.g. M 82, Thuma et al. l2TjtjU|) . There the line 
contribution to the measured continuum flux ranges from 
8 % in NGC 5907 to 35 % in the centre of NGC 3079, and 
it seems to be correlated with the activity status of the 
galaxy. This would point to NGC 4631 being a rather inac- 
tive galaxy. However, a small fraction of line emission can 
also point to a higher value for the dust emission per gas 
mass, i.e. to a higher absorption cross section (see further 
below). 

The two measured flux density values at A A 0.87 and 
1.23mm, together with the IRAS data (Young et al. lT^ll . 
still leave a large gap in the spectral coverage between the 
FIR and the sub-mm. In order to fill this gap, we take the 
ISO and SCUBA data published by Bendo et al. 
2003) and estimate flux density values for wavelengths of 
180 fim and 450 /im. 

Concerning the SCUBA data, these authors give a flux 
density value only for the inner 135" of the galaxy. In or- 
der to estimate the total flux density for NGC 4631, we 
have to scale this by a factor S , tot/<5'i35" , describing the 
ratio of the total flux density to the flux density in the 
central 135" of the galaxy. From our maps, we measured 
this ratio to 2.27 and 2.34 at a wavelength of 0.87mm 
and 1.23 mm, respectively (after subtraction of the non- 
dust contributions, which are negligible at A450/im). In 
the FIR, the ratio of IRAS flux densities (corresponding to 
the total emission) to the ISO values (corresponding to the 
inner 135") is about 1.6 at A60/xm and 1.5 at AlOO/xm, 
where the latter value is less certain due to the wider 
PSF of ISOPHOT at this wavelength. Considering that 
the warmer dust (i.e. at shorter wavelengths) is probably 
more concentrated to the inner disk than the colder dust, 
we assume a factor of <Stot /5i35" = 2.0 at A 450 /im. With 
this value we get a total flux density of S^o^m = 36 ±9 Jy. 

In order to estimate the total flux density at A 180 /im, 
wc take into account that the ISOPHOT C200 detector 
array covers only the central 180" x 180" of the galaxy, 
and therefore also misses a significant fraction of the 
disk emission. Again we can estimate a scaling factor 
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Fig. 4. The FIR-to-mm spectrum of NGC 4631. The FIR 
data points up to A 100 /im are taken from Young et al. 
(1989), those at 180 /im and 450 /im are estimated from 
the values in Bendo et al. (2QJE3, and the 0.87 mm and 
1.23 mm flux densities are from this paper and corrected 
for line contributions and synchrotron and free-free emis- 
sion. The solid line shows a fit for a "standard" dust model 
(two modified Planck curves with = 2). 



from our maps, and find Stot/Siso" ~ 1.8 for A180/im. 
Interestingly, we find at AA 0.87 mm and 1.23 mm that the 
flux density within the central 180" is not much larger 
than that within the central 135". This can be understood 
by the fact that the smaller 135" area covers already the 
central "dust ring", as described in the previous section. 
This result differs from the values given by Bendo et al. 
(2003J), who found a much larger difference between the 
two flux density values for the inner 135" and the inner 
180". This inconsistency can probably be explained by the 
fact that these authors used a SCUBA sub-mm map for a 
deconvolution analysis, although this map covers less then 
3' of the inner disk, and is therefore not well suited for this 
purpose. 

In addition, the spatial resolution of ISO at this wave- 
length is not much smaller than the detector array. With 
the knowledge of the PSF of this detector one can deter- 
mine which fraction of the total emission of a source would 
be detected by the C200. For a somewhat extended source 
like NGC 4631 we estimate a value of 82%. Taking these 
two corrections into account, we get a total flux density 
at this wavelength of Siso^m = 205 ± 35 Jy. 
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4.2. Dust temperatures and composition 

The interstellar dust in galaxies is a mixture of many 
components at several temperatures. Although the exact 
composition and temperature distribution depends on the 
local interstellar radiation field and local processes im- 
portant for dust grain processing, the total FIR to mm- 
spectrum of a galaxy can usually be fitted well by 1 - 3 
components with a modified black-body spectrum of the 
form 



S\ oc erf B\(T d ) , 



(5) 



with the dust absorption cross section per hydrogen atom 
ct^ oc A"' 3 . The value of /3 depends on the dust com- 
position and is still a matter of debate; theoretical ar- 
guments suggest values between 1 and 2. Realistic dust 
models, which are able to explain the observed dust emis- 
sion and extinction, suggest (3 ~ 2 for big grains which 
radiate in thermal equilibrium with the ambient temper- 
ature (Andriesse 1974: Draine & Lee 1984). In fact the 
FIR-to-mm spectra of many external galaxies have been 
successfully fitted with one or two components of dust 
with (5 = 2 (e.g. Chini et al. IT5351 Neininger et al. 11391)1 
Dumke et al. ITMTI Braine et al. [TWH Alton et al. IHMl 
Stevens & Oea.r IMTO . 

We fitted the FIR-to-mm spectrum of NGC 4631 with 
a dust model consisting of two components, both with 
(3 = 2. We used data points from Young et al. (1989J in 
the FIR (25 - 100 /xm), those estimated above at 180 /mi 
and 450 /jm from Bendo et al. (2003|), and our new values 
at 0.87 mm and 1.23 mm. The result is shown in Fig. 
where the solid curve represents a two-component modi- 
fied Planck spectrum with temperatures of 50 K and 18 K. 
This temperatures agree well with previous results (see 
references above) , where the cold dust component was usu- 
ally found to have a temperature of 15 - 20 K. 

We can also estimate dust temperatures for the inner 
and outer disk separately when we use the flux density 
values at AA60/im and 100 fim from Bendo et al. (2003) 
for the inner 135", and restrict the integration in our maps 
to the same area. A fit to the data leads to somewhat 
higher dust temperatures of 51 K and 20 K for the warm 
and cold component, respectively. After subtracting the 
flux density values for the inner disk from the total flux 
densities, we also estimate dust temperatures for the outer 
disk and find 49 K for the warm component and 16 K for 
the cold component. Please note, however, that this two- 
component model is still a model, and the main conclusion 
from these fits is that the dust temperature decreases from 
the centre of NGC 4631 to the outer disk. 

As it can be seen in Fig.0] this standard model seems 
to be not sufficient to explain the measured flux densities 
at both 0.87 and 1.23 mm. Instead we detect an excess 
in the mm-range (or, alternatively, a lack of emission in 
the sub- mm): The ratio of the two dust-only flux densi- 
ties (i.e. after subtraction of line and other contributions) 
is S , dust(0.87mm)/S , dust (1.23mm) = 1.67±0.35, while we 
would expect a ratio of 4 in case of a modified Planck 



spectrum with a dust absorption coefficient oc A~ 2 . This 
excess of dust emission at 1.23 mm leads to a flattening of 
the sub- mm/mm spectrum of NGC 4631, which is difficult 
to explain with the thermal radiation of big dust grains at 
moderately warm temperatures alone. However, for prac- 
tical reasons, we continue here with the analysis of the 
data by assuming that the above model is reasonable, and 
discuss this mm excess further below. 

4.3. Absorption cross sections 

An averaged value of the dust absorption cross section 
in NGC 4631 can be calculated from the measured flux 
density and the hydrogen mass through (e.g. Hildebrand 
1983) 



Sx m H D 2 
B x (T d )M H 



(6) 



Here mjj is the mass of a hydrogen atom, D the distance 
to the galaxy, and Mh = Mhi + Mh 2 the hydrogen mass. 
From observations of the atomic hydrogen in NGC 4631, 
Rand (1994) finds a Hi mass of the disk (excluding the 
Helium content) of Mm = 5 x 10 9 M©. We don't take into 
account here the mass of the various Hi spurs detected 
by Rand (|1994p . which would add another 2.2 x 10 9 M Q . 
The molecular mass found by Golla & Wielebinski (1994) 
is Mh 2 = 1.2 x 10 9 M©, based on the galactic CO-H 2 con- 
version factor of X = 2.3 x 10 20 cm _2 /(Kkms~ 1 ) (Strong 
et al. I1988|) . Although we know that X may be variable 
and not necessarily valid for NGC 4631, the fraction of 
the molecular mass is rather small in this galaxy, thus we 
use M H = 6.2 x 10 9 M Q . When we take this value for 
the hydrogen mass, assume a dust temperature of 18 K 
as determined above (note that at our observing wave- 
lengths around 1 mm more than 99 % of the dust emission 
comes from the cold component; see Fig.0J, and use Eq. 
© to calculate the absorption cross section for the whole 
galaxy, we get values of cr^ 8 7 mm = 6.2 ± 1.0 x 10~ 27 cm 2 
and cr^mm = 6.4 ±0.8 x 10~ 27 cm 2 for both wavelengths. 

In order to interpret these values, we use a formulation 
following Mezger et al. I|1990(l . 



(7) 



C = 7 ■ 10 -27 cm 2 is the dust absorption cross section at 
A 1 mm following from the theoretical curves of Draine & 
Lee (|1984fl . A mm the wavelength in mm, Z the metallicity 
and b an empirically determined factor which accounts for 
the differences between Drain & Lee's grain mixtures and 
real grains. 6=1 applies to dust in the diffuse interstel- 
lar Hi gas, and since most gas in NGC 4631 is in atomic 
form, we will also use b = 1 here. For solar metallicity we 
would thus expect values of 9.2 and 4.6xl0 _27 cm 2 for 
0.87 and 1.23 mm, respectively. At the latter wavelength, 
such a value has been measured for several nearby galaxies 
(e.g. Neininger et al. lTM)! Dumke et al UMTjl . Kriigel & 
Chini (1994) give a "standard" value which corresponds 
5.0 x 10 _27 cm 2 based on typical abundances 
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of silicate and graphite in the ISM, which is in very good 
agreement with the values above and has been used suc- 
cessfully for fitting IR spectra over the whole spectrum 
from 1 to 1300 /im. 

To compare our results with these values, we have to 
take the metallicity into account: Otte et al. I|2002[) mea- 
sured the Nitrogen abundance in NGC 4631 and found an 
average metallicity of Z/Zq ~ 0.5. This is also in agree- 
ment with the results from Vila-Costas & Edmunds (1992) 
who determined metallicities and their gradients in Sd 
galaxies. With Z = 0.5 Zq, the observed absorption cross 
section a is a factor 1.5 higher than the predicted value for 
A 0.87 mm, and even a factor 3 higher than the predicted 
value for A 1.23 mm. 

The dust absorption cross section can also be expressed 
by (Hildebrand 



ax = k\ toh 



Md 



(8) 



with the dust absorption coefficient n\ (which involves 
grain properties like size, density, and emissivity) and 
the dust-to-gas mass ratio Md/M g . The high values for 
a can thus be explained either with a high dust-to-gas 
ratio (which may be in contradiction to the low metal- 
licity) or an increased absorption coefficient in the sub- 
mm/mm regime, and therefore unusual optical properties 
of the dust grains in NGC 4631. 

From Eq. JHJ it is obvious that these high values for 
a are due to rather high flux densities. But when com- 
paring the numbers for both wavelengths, the high value 
at A 1.23 mm can also be directly related to the unusual 
ratio between the two flux densities at our observing wave- 
lengths. On the other hand, as shown in Fig. we also 
see that the intensity ratio between these two wavelengths 
varies strongly along the galaxy's major axis. Therefore 
we also calculate the absorption cross section locally for 
selected areas along the major axis of NGC 4631, in order 
to trace local variations of the unusual dust properties 
found in the total spectrum. 

The flux density per beam emitted by a cloud of gas 
and dust is (see Mezger et al. [1990) 



5 x = n b e am -BA(T d )(l-e-^) 



(9) 



with t\ = o-^Nh, where is the dust absorption cross 
section per hydrogen atom and iVu the beam-averaged 
hydrogen column density. For rj < 1 we can transform 
Eq. © and express the cross section by 



A 2 



beam 



2kT d 



_1) jSa 



(10) 



with x = hc/XkT d and N H = iV(Hi) ± 2 7V(H 2 ). 

We can use Eq. (fT~0f> to calculate the absorption cross 
section in various parts of NGC 4631. For this calculation 
we assume the temperature of the cold dust to decrease 
linearly from 20 K in the centre to 16 K in the outer parts, 
as found in the previous subsection. The result is shown in 
the upper panel of Fig. [51 which also shows (in the lower 




400 



x offset [arcseo] 



Fig. 5. Distribution of the intensity ratio of the two ob- 
serving wavelengths (lower panel) and the absorption cross 
sections for a standard dust model (upper panel) along the 
major axis of NGC4631. 

Table 3. Dust emission intensity ratios and absorption 
cross section in various regions along the major axis of 
NGC 4631. 



Position 


^0.87mm/^1.23mm 


H 

°0.87mm 

[xl(T 


H 

"1.23mm 

27 cm 2 ] 


x > 4' (east) 


2.8 ± 0.6 


15 ±4 


8± 2 


4' > x > 1'.5 


1.3 ±0.2 


7± 2 


9± 1 


1'.5 > x > -lf5 


1.7 ±0.1 


7± 1 


7± 1 


-1'.5 > x > -4' 


1.4 ±0.1 


9±2 


12 ± 1 


x < —4' (west) 


1.7 ±0.2 


19 ±3 


18 ±2 



panel) the variation of the intensity ratio Io.s7mm/ h. 23mm 
along the major axis. 

Before we discuss these data, we average the results 
in various parts of the major axis in order to increase the 
S/N: in the centre (strong CO, |a;| < 1?5), in the inner 
disk (weak CO, 1'.5 < \x\ < 4'), and in the outer disk (no 
detected CO, \x\ > 4'). The obtained numbers for the dust 
intensity ratios and absorption cross sections are given in 
Table El 

From these numbers, we can roughly distinguish three 
different situations. In the centre and inner disk, the ab- 
sorption cross section at A 0.87 mm is only slightly higher 
than expected, whereas it is much higher at A 1.23 mm. 
Also the intensity ratio is very low, especially in the areas 
of weak CO emission, where 01.23mm is even higher than 
Co. 87mm- In the outer parts however, where the optical 
image of NGC 4631 shows strong disturbances, the situ- 
ation has changed. At the eastern edge of the disk, also 
Vo.srmm ^ s strongly enhanced, while the intensity ratio be- 
tween both wavelengths is closer to the expected values. 
At the western edge, on the other hand, the intensity ratio 
is similar to the central area, but at both wavelengths the 
absorption cross section is increased by another factor 2 
compared to the central area and inner disk. 

To summarize the results of this subsection, we find 
that the dust absorption cross section in the sub- mm/mm 
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range is higher than what is predicted, and this result 
is especially significant in the outer disk (\x\ > 4') of 
NGC4631, where no molecular gas is detected. These ar- 
eas are also strongly affected by the gravitational inter- 
action of NGC 4631 with its neighbours. The optical im- 
age reveals strong disturbances of the disk, and also the 
footpoints of two Hi spurs detected by Rand (1994J are 
located in these areas. A significant amount of molecular 
gas in these regions without corresponding CO emission 
could be responsible for the higher values of the absorp- 
tion cross section at A 0.87 mm, since this would result in 
an underestimate of Nr and therefore a too high value of 
er H (see Eq. 11100 ). This possibility is supported by results 
which suggest a higher value for the CO-H2 conversion fac- 
tor in the outer disks of galaxies and metal-poor environ- 
ments - note that the metallicity in NGC 4631 decreases 
with increasing galactocentric radius (Otte et al. I2002[l . 
Nevertheless this possibility cannot explain why cri.23mm 
is of the same order or even larger than co.87mm, or the ob- 
served excess of the dust emission at A 1.23 mm. Another, 
more speculative scenario, would involve a different way 
of dust processing in those areas of NGC 4631 which are 
strongly affected by the interaction. 

Therefore a reasonable dust model, whose properties 
can account for the mm excess detected above, must 
also be able to explain the enhanced values of a at sub- 
mm/mm wavelengths. 

4.4. What is the origin of the mm excess? 

As mentioned above, the sub-mm/mm spectrum of 
NGC 4631 is too flat to be satisfactorily explained with a 
"standard" two-component dust model. However, several 
other dust models have been proposed in the past, includ- 
ing stochastic heating of very small particles, big grains at 
very low temperatures (<C 10 K), or grains with different 
optical properties, like fractal, fluffy, or ice-coated grains. 
A good overview over several models is given by Reach et 
al. I|1995(l in an attempt to interpret COBE observations 
of our Galaxy. In the following we summarize the main 
points of these models in view of their application to the 
case of NGC 4631. 

4.4.1. Very cold big grains 

In most external galaxies which have been observed in the 
mm and sub-mm continuum the cold dust component has 
a temperature of 15 - 20 K. However, Siebenmorgen et al. 
(1999 ) have detected very cold dust in a sample of inac- 
tive spiral galaxies with an average temperature of ~ 13 K. 
These authors also discuss the possible existence of even 
colder dust and conclude that within the optical disk of 
a galaxy the dust cannot be colder than about 6 K. If we 
want to explain the mm-excess of NGC 4631 with such a 
component of very cold grains, the most reasonable so- 
lution yields a temperature of T vc = 4K for this tem- 
perature component, and the upper limit we find is 6K. 



Such a component has its maximum close to A ~ 1 mm 
and could therefore account for the very flat sub-mm/mm 
part of the flux density spectrum. The fit to the data with 
temperatures of 50 K, 20 K, and 4K is shown in Fig. 

Such very cold dust could in principle exist in the form 
of self-shielded grains in very optically thick clouds with 
no intrinsic sources. In order to produce the measured flux 
density at A 1.23 mm, the total dust mass needed is of the 
order 3 — 6 x 10 8 M Q (assuming an absorption coefficient 
as in the solar neighbourhood), most of which is at a tem- 
perature below 6 K. With the gas mass estimated from the 
CO and Hi data (8.5 x 10 9 M Q including the Helium con- 
tent) this leads to a gas-to-dust ratio of 15 — 30, far below 
values in other galaxies or the Milky Way. While this low 
ratio could explain the high value we found for the dust 
absorption cross section a (see Eq. ©), it is in contradic- 
tion to the low metallicities in NGC 4631 (Z/Z Q ~ 0.5 on 
average) which were found recently by Otte et al. ( 2002 ) . 

On the other hand, Kriigel & Siebenmorgen (1994) 
have shown that in cold dense clouds the dust grains may 
coagulate to ice-coated very big grains (up to a radius 
a ~ 100 /mi), for which the absorption coefficient can be 
enhanced by a factor of eight relative to the diffuse ISM. 
If a large amount of cold dust existed in this form, the 
dust emission could be explained with gas-to-dust ratios 
similar to values in the local ISM. 

However, in order to keep a large amount of dust at 
temperatures below 6K, compared to the cold compo- 
nent with 15 - 20 K, the interstellar radiation field must 
be significantly attenuated. The necessary extinction is 
of the order Ay = 15 — 50 (depending on grain type), 
corresponding to a gas column density of iVn ~ 10 22 or 
more (see Reach et al. 119951 and references therein). Fig. El 
shows that values of Nh ~ 10 22 for the total hydrogen col- 
umn density are reached only in the inner 3' of NGC 4631, 
around the maxima of the CO distribution (assuming a 
standard CO-H 2 conversion factor). While we would ex- 
pect radiation-shielded cold dust to exist at places of dense 
molecular gas, these regions are also places of strong star 
formation. In addition, the mm excess is present over the 
whole disk of NGC 4631, and especially at radii of lf5 — 3', 
where the ratio Jo. 87mm /-^l. 23mm is smallest (see Fig. |SJ). 
These areas are located well beyond the molecular ring. 

On the other hand, Galliano et al. (2003] consistently 
explained the mm-excess in NGC 1569 with very cold 
grains deeply embedded in clumps of big grains and pri- 
marily heated by the FIR emission of the latter. This 
would require a very clumpy medium and small filling fac- 
tors of the cold gas and dust. 

4.4.2. Very small grains 

A different model which was successfully applied to the 
absorption and emission data from our galaxy over a 
large wavelength range was presented by Desert et al. 
(1990). This model included (besides big grains, which 
radiate mainly at A > 70 /im) a component of amorphous, 
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Fig. 6. Two alternative possibilities to fit the FIR-to-mm 
spectrum of NGC 4631. The left figure shows a possible in- 
terpretation including a component of very cold dust. The 
right figure shows a model fit with an increased abundance 
of very small grains. Both possibilities are discussed in the 
text. 



very small grains (VSGs), which are the main contrib- 
utor to the IR emission between 10 and 70 /im, and of 
Polycyclic Aromatic Hydrocarbons (PAHs), which radi- 
ate shortwards of 10 /im. The VSGs, which are about 1 - 
20 nm in size, are so small that they are heated to non- 
equilibrium temperatures by absorption of a single pho- 
ton and thus show a fluctuating temperature distribution. 
They also have a much broader FIR emission spectrum 
compared to big grains, and an absorption coefficient pro- 
portional to A -1 , i.e. (3—1 can be assumed for this com- 
ponent (Seki & Yamamoto 1980). 

Besides our own galaxy, also the FIR-to-mm spectrum 
of the dwarf galaxy NGC 1569 was succesfully fitted with 
this dust model (Lisenfeld et al. I2002|l , although with a 
stronger ambient radiation field and a different grain com- 
position (increased abundance of VSGs and absence of 
PAHs) than the solar neighbourhood. 

We tried to fit the FIR to mm spectrum of NGC 4631 
with the dust model of Desert et al. I|1990(l . the result is 
shown in Fig. [BJ For this fit we neglected the PAH compo- 
nent which radiates at wavelengths smaller than 25 /im. In 
order to fit the data, we have to assume an interstellar ra- 
diation field and dust composition different from the solar 
neighbourhood. From the FUV data published by Smith et 
al. H2001[> . we estimate an average ISRF within NGC 4631 
at A 1000 A of about 6 x 10" 2 erg cm 2 s /im , which 
is four times the local ISRF as given in Mezger et al. 
(1982). Compared to big grains, which have an average 
temperature of ~ 22 K in this model, VSGs are over- 
abundant by a factor of 4 in order to account for the 
small S'(0.87mm)/S'(1.23mm) ratio with the shallower 
VSG spectrum in the FIR. In addition, we assume that 
the VSGs have sizes up to 20 nm, compared to 15 nm as 
in the Galaxy. This latter assumption does not necessarily 
change the physical properties of the VSGs, like surface- 
to- volume ratio and heat capacity (Kamijo et al. Jl975 
Stephens & Russell EMi . 

The application of this model to NGC 4631 has its 
drawbacks as well. The estimated dust mass for the big 



grain component, assuming standard optical properties, is 
1.1 x 10 7 M Q , while the VSGs add up to 0.3 x 1O 7 M . The 
resulting gas-to-dust ratio of ~ 600 is similar to values in 
extremely metal-poor environments, and therefore some- 
what too high for the metallicities measured in NGC 4631. 
Another critical point is the ISRF of 4 times solar which 
is required for a reasonable fit. While this value seems ap- 
propriate in the inner part of NGC 4631, it is unclear if 
and how such a radiation field can be maintained in the 
outer disk, where the measured absorption cross sections 
are highest. Besides these physical concerns, the model fit 
to the data resulted in a x 2 much worse than the model 
including a component of very cold grains. 

4.4.3. Grains with unusual optical properties 

Fractal grains have an enhanced efficiency at sub-mm 
and mm-wavelengths compared to the optical (e.g. Wright 
1993). Thus a population of fractal grains may exist which 
show rapid temperature fluctuations (of a few K), but 
spend most of the time at very low temperatures due 
to high FIR-to-mm emissivity. While the exact proper- 
ties of this type of grains are not yet well established, 
it is known that a much smaller amount of these grains 
is needed (compared to compact spherical grains) in or- 
der to produce a comparable amount of emission, hereby 
explaining the large values for a found in the previous 
subsection. And while the temperature spectrum of these 
grains may well account for the measured mm flux density, 
their mass plays only a minor role in the total dust mass 
of the galaxy. Thus the gas-to-dust ratio could have a rea- 
sonable value and agree with the low metallicity inferred 
from the Nitrogen abundance. 

If such a component exists, we have to ask why it is 
detected only in a few objects. It may be responsible for 
part of the mm radiation of the Milky Way (Reach et al. 
1995) and the dwarf starburst galaxy NGC 1569 (Galliano 
et al.EnnHI Lisenfeld et al.EEEI), although the latter au- 
thors explained this emission with an overabundance of 
VSGs. Many other galaxies do not show any excess in the 
mm range, but why should this fractal grain population 
exist in these apparently very different objects, but not 
in many other galaxies of various types? We can argue 
that sub-mm continuum observations have just reached 
the required sensitivity to detect this effect, and for many 
objects only data in the mm or in the sub-mm regime ex- 
ist, while our results show that it is necessary to observe 
at both wavelengths to investigate the dust properties of 
these objects. 

After all we should note that in principle strong vari- 
ations in strength and colour of the interstellar radia- 
tion field may also cause a broadening of the dust spec- 
trum, hereby potentially flattening the spectrum in the 
sub- mm/mm range. However, shortwards of A 0.87 mm the 
spectrum is not flattened and consistent with a dust ab- 
sorption coefficient cx A~ 2 . Furthermore, the discovered 
excess in the mm is quite significant, and with the com- 
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plexity of the ISM also in normal galaxies the question 
arises why the temperature distribution of dust should be 
so much narrower in most other galaxies observed up to 
date, despite their varying classification and star forma- 
tion properties. 

5. Summary and outlook 

We have observed the nearby interacting galaxy NGC 4631 
in the radio continuum emission at AA 0.87 mm and 
1.23 mm. The emission is concentrated on the galactic 
plane, although some halo emission could also be detected, 
especially at 1.23 mm where intergalactic cold dust seems 
to follow the previously detected H I spurs. The two max- 
ima visible in both maps, located in the plane symmetri- 
cally around the nucleus at a radius of about 30", resemble 
a ring-like structure, similar to the distribution of the CO 
molecular line emission. In the outer disk, where no CO 
is detected, the mm/sub-mm emission follows the Hi dis- 
tribution, similar to other edge-on galaxies. 

After subtracting non-dust contributions (line, syn- 
chrotron, and free-free emission) from the bolometer flux, 
the remaining flux densities of So. 87mm = 3.41 ± 0.58 mJy 
and Si. 23mm = 2.04 ± 0.23 mJy can be attributed to the 
thermal emission of dust. In order to estimate the tem- 
perature of this dust, we fitted a two-component modified 
Planck spectrum to the FIR-to-mm data and found a tem- 
perature of 50 K for the warm component and of 18 K for 
the cold component, which is responsible for most of the 
emission at A > 200 /zm. 

This two-component model suffers from the fact that 
the observed dust spectrum is too flat in the sub-mm/mm 
range. Furthermore it leads to absorption cross sections 
too high compared with theoretical expectations or previ- 
ous results on other galaxies, especially in the outer part 
of the disk which are disturbed because of gravitational 
interaction. 

We suggest several possibilitcs to resolve these incon- 
sistencies, the most likely being grains with unusual op- 
tical properties which can account for the measured mm- 
excess as well as the high absorption cross sections. 

These observations have shown that results of mm and 
sub-mm observations have to be combined to investigate 
the cold dust component in external galaxies. 
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